Measurements of G^/G^ from the 2 H(e,e'n) 1 H Reaction to Q 2 = 1.45 (GeV/c) 2 



R. Madey, 1 ' 2 A. Yu. Semenov, 1 S. Taylor, 3 B. Plaster, 3 A. Aghalaryan, 4 E. Crouse, 5 G. MacLachlan, 6 S. Tajima, 7 
W. Tireman, 1 Chenyu Yan, 1 A. Ahmidouch, 8 B. D. Anderson, 1 H. Arenhovel, 9 R. Asaturyan, 4 O. Baker, 10 
A. R. Baldwin, 1 D. Barkhuff, 3 'E H. Breuer, 11 R. Carlini, 2 E. Christy, 10 S. Churchwell, 7 'Q L. Cole, 10 
S. Danagoulian, 2 ' 8 D. Day, 12 T. Eden, 1 ' lo <0 M. Elaasar, 13 R. Ent, 2 M. Farkhondeh, 3 H. Fenker, 2 J. M. Finn, 5 

L. Gan, 10 K. Garrow, 2 P. Gueye, 10 C. R. Howell, 7 B. Hu, 10 M. K. Jones, 2 J. J. Kelly, 11 C. Keppel, 10 
M. Khandaker, 14 W.-Y. Kim, 15 S. Kowalski, 3 A. Lai, 1 A. Lung, 2 D. Mack, 2 D. M. Manley, 1 P. Markowitz, 16 

J. Mitchell, 2 H. Mkrtchyan, 4 A. K. Opper, 6 C. Perdrisat, 5 V. Punjabi, 14 B. Raue, 16 T. Reichelt, 17 
J. Reinhold, 16 J. Roche, 5 Y. Sato, 10 N. Savvinov, 11 I. A. Semenova, 1 W. Seo, 15 N. Simicevic, 18 G. Smith, 2 
S. Stepanyan, 4 * 15 V. Tadevosyan, 4 L. Tang, 10 P. E. Ulmer, 19 W. Vulcan, 2 J. W. Watson, 1 S. Wells, 18 
F. Wesselmann, 12 S. Wood, 2 Chen Yan, 2 S. Yang, 15 L. Yuan, 10 W.-M. Zhang, 1 H. Zhu, 12 and X. Zhu 10 

(The Jefferson Laboratory E93-038 Collaboration) 

1 Kent State University, Kent, Ohio 44^4^ 
2 Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 
J Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
* Yerevan Physics Institute, Yerevan 375036, Armenia 
The College of William and Mary, Williamsburg, Virginia 23187 

6 Ohio University, Athens, Ohio 45701 
7 Duke University and TUNL, Durham, North Carolina 27708 
8 North Carolina A&T State University, Greensboro, North Carolina 27411 
Johannes Gutenberg- Universitat, D-55099 Mainz, Germany 
10 Hampton University, Hampton, Virginia, 23668 

11 University of Maryland, College Park, Maryland 20742 

12 University of Virginia, Charlottesville, Virginia 22904 

1S Southern University at New Orleans, New Orleans, Louisiana 70126 
'^Norfolk State University, Norfolk, Virginia 23504 
15 Kyungpook National University, Taegu 702-701, Korea 
16 Florida International University, Miami, Florida 33199 
17 Rheinische Friedrich-Wilhelms-Universitat, D-53115 Bonn, Germany 
1 Louisiana Tech University, Ruston, Louisiana 71272 
19 Old Dominion University, Norfolk, Virginia 23529 
(Dated: February 8, 2008) 

We report new measurements of the ratio of the electric form factor to the magnetic form factor of 
the neutron, G^/GJf, obtained via recoil polarimetry from the quasielastic 2 H(e, e'n) 1 !! reaction at 
Q 2 values of 0.45, 1.13, and 1.45 (GeV/c) 2 with relative statistical uncertainties of 7.6 and 8.4% at 
the two higher Q 2 points, which were not reached previously via polarization measurements. Scale 
and systematic uncertainties are small. 

PACS numbers: 14.20.Dh, 13.40.Gp, 25.30.Bf, 24.70,+s 



The nucleon elastic electromagnetic form factors are 
fundamental quantities needed for an understanding of 
nucleon and nuclear structure. The evolution of the elec- 
tric and magnetic form factors with Q 2 , the square of the 
four-momentum transfer, is related to the charge and cur- 
rent distributions within the nucleon. Precision measure- 
ments of the electromagnetic form factors are important 
for tests of non-perturbative quantum chromodynamics 
(QCD) either on the lattice or in models. Measurements 
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of the neutron electric form factor, G^., have been im- 
peded by the lack of a free neutron target and the small 
value of G\ relative to the neutron magnetic form factor, 
G 7 ^; however, with the advent of high duty- factor polar- 
ized electron beam facilities, experiments employing re- 
coil polarimeters 0,0], polarized 3 He targets jUGIB, an< ^ 
polarized deuterium targets 0, 0] have yielded the first 
precision measurements of G^>. These polarization mea- 
surements of Gg are limited to Q 2 < 0.67 (GeV/c) 2 and 
are, within errors, consistent with the Galster parame- 
terization |8(. In addition to polarization measurements, 
G^ has been extracted from a theoretical analysis of the 
deuteron quadrupole form factor @ for Q 2 values up to 
1.6 (GeV/c) 2 . 

In the plane-wave approximation, the recoil polar- 
ization produced by a longitudinally polarized electron 
beam in quasielastic electron-neutron scattering is re- 



2 



Top Rear Array 




Rear Veto/Tagger 

/Front Array 
j 



Bottom Rear Array 

Front Veto/Tagger 





Lead Curtain C nar ybdis 

Target LD2, LH2 

FIG. 1: A schematic diagram of the polarimeter. 

stricted to the scattering plane [T(llrTI |: The longitudinal 
component, Pjj , and the transverse (sideways) compo- 
nent, Ps>, are parallel and perpendicular, respectively, 
to the recoil neutron's momentum vector. In terms of 
and G^, Ps> and Py can be written as 



Ps'/Pl 
Pl'/Pl 



-KsG n E G n M i 'I , 
K L {G\tflh , 



(1) 
(2) 



where Pl is the electron beam polarization, Iq=(G%) 2 + 
Ko(Gj^) 2 , and K$, Kl, and Kq are kinematic functions 
of the electron scattering angle, 6 e , and Q 2 . Measure- 
ments of Ps' and Py via a secondary analyzing reaction 
permit an extraction of the ratio of Gg to G^; a sig- 
nificant advantage of this technique is that Pl and the 
analyzing power of the secondary reaction cancel in the 
polarization ratio Ps'/Pl'- Also, for quasifree emission, 
Arenhovel ^3 demonstrated that P$> and Py are insen- 
sitive to final state interactions (FSI), meson exchange 
currents (MEC), isobar configurations (IC), and to the- 
oretical models of deuteron structure. In this Letter, we 
report new measurements of G^/G^- obtained via re- 
coil polarimetry from the quasielastic 2 H(e, e'n) 1 !! reac- 
tion at three central Q 2 values of 0.45, 1.15, and 1.47 
(GeV/c) 2 . 

Our measurements were carried out in Hall C of the 
Thomas Jefferson National Accelerator Facility. The ex- 
perimental arrangement with an isometric view of our 
polarimeter is shown in Fig.^ A beam of longitudinally 
polarized electrons (with a typical polarization of 80%) 
scattered quasielastically from a neutron in a 15-cm liq- 
uid deuterium target. A scattered electron was detected 
in the High Momentum Spectrometer (HMS) in coinci- 
dence with the recoil neutron. The neutron polarime- 
ter (NPOL) was used to measure the up-down scattering 
asymmetry from the transverse component of the recoil 
neutron polarization presented to the polarimeter. To 
permit measurements of the up-down scattering asymme- 
try from different combinations of Ps> and Py , a dipolc 
magnet (Charybdis) located in front of the polarimeter 
precessed the recoil neutron's polarization vector through 
an angle \- 
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FIG. 2: Typical time-of-flight spectra for Q 2 = 1.15 
(GeV/e) 2 . The dark-shaded regions indicate the selected por- 
tions of the spectra. 



The polarimeter consisted of a total of 44 plastic scin- 
tillation detectors. To achieve luminosities of ~ 3xl0 38 
cm~ 2 s _1 , the front array was segmented into 20 detec- 
tors [100 cmxlO cmxlO cm]. Top and bottom rear ar- 
rays were shielded from the direct path of particles from 
the target. Each rear array consisted of 6 "20- in" de- 
tectors [101.6 cmx50.8 cmx 10.16 cm] and 6 "10-in" de- 
tectors [101.6 cmx25.4 cmxl0.16 cm]. A double layer 
of "veto/tagger" detectors (each 0.64-cm thick) directly 
ahead of and behind the front array identified incom- 
ing and scattered charged particles. A 10-cm lead cur- 
tain attenuated the flux of electromagnetic radiation and 
charged particles incident on the polarimeter. The flight 
path from the center of the target to the center of the 
front array was 7.0 m, and the mean flight path from the 
front array to the rear array was 2.5 m. 

For a fixed neutron scattering angle of 46.0°, central 
Q 2 values of 0.45 and 1.47 (GeV/c) 2 were associated with 
beam energies of 0.884 and 3.40 GeV, respectively, and 
electron scattering angles of 52.7° and 23.6°, respectively. 
The measurement conducted at a central Q 2 value of 
1.15 (GeV/c) 2 was associated with two beam energies 
of 2.33 and 2.42 GeV and electron scattering angles of 
30.8° and 30.1°, respectively. We conducted asymme- 
try measurements with the polarization vector precessed 
through x = ±40° at each of our Q 2 points; in addition, 
at Q 2 = 1.15 and 1.47 (GeV/c) 2 , we conducted asymme- 
try measurements with the polarization vector precessed 
through x — 0°> ±90°. The acceptance- averaged values 
of Q 2 are: (Q 2 ) = 0.45, 1.13, and 1.45 (GeV/c) 2 . 

Typical time-of-flight spectra are shown in Fig. [21 The 
left panel is an HMS-NPOL coincidence time-of-flight 
spectrum. We compared the measured time-of-flight, 
cTOF, with the time-of-flight calculated from electron 
kinematics and offsets determined by a calibration pro- 
cedure; the result is centered on zero with a FWHM of 
approximately 1.5 ns. The right panel is the time-of-flight 
spectrum between a neutron event in the front array and 
an event in the top or bottom rear array. We compared 
this measured time-of-flight, ATOF, with the time-of- 
flight calculated for elastic np scattering. Variations with 
respect to a nominal 2.5 m flight path were compensated. 
The tail on the slow side is due to Fermi motion in car- 
bon and nuclear reactions, and the secondary peak at 
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— — 2.5 ns is the result of ir° production in the front array. 
To extract the physical scattering asymmetry, we calcu- 
lated the cross ratio, r, which is defined to be the ratio 
of two geometric means, (A+A^j) 1 / 2 and (N^N+) 1 ^ 2 7 
where Njj(Np) is the yield in the ATOF peak for neu- 
trons scattered up (down) when the beam helicity was 
positive(negative); the yields, corrected for background, 
were obtained by peak fitting. The physical scattering 
asymmetry is then given by (r — l)/(r + 1). The merit 
of the cross ratio technique |l3j is that the neutron po- 
larimeter results are independent of the luminosities for 
positive and negative helicities, and the efficiencies and 
acceptances of the top and bottom halves of the polarime- 
ter. Beam charge asymmetries (of typically 0.1%) and 
detector threshold differences cancel in the cross ratio. 

To account for the finite experimental acceptance and 
nuclear physics effects such as FSI, MEC, and IC, we av- 
eraged Arenhovel's theoretical 2 H(e,e'n) 1 H calculations 
|14| over the experimental acceptance. These calcula- 
tions include leading-order relativistic contributions to a 
non-relativistic model of the deuteron as an n-p system, 
employ the Bonn i?-space NN potential ^| for the in- 
clusion of FSI, and include MEC and IC. Other realistic 
potentials (e.g., the Argonne V18 0|) give essentially the 
same results. The theoretical values of the recoil polariza- 
tion were calculated over a kinematic grid; during the av- 
eraging procedure, the recoil polarization was computed 
via multidimensional interpolation between the grid ele- 
ments. 

To average these theoretical calculations over the ex- 
perimental acceptance, we prepared two independent 
simulation programs. First, we developed the GENGEN 
Monte Carlo simulation program, which includes an 
event generator and detailed models of the electron spec- 
trometer and the neutron polarimetcr. GENGEN repro- 
duces experimental kinematic distributions and models 
the response of the polarimeter. Second, we developed 
a program that used the kinematics of the reconstructed 
quasielastic events from the experimental data to com- 
pute the recoil polarization for each event used in the 
data analysis; the advantage of this method is that it 
does not require a model of the experimental acceptance. 

For the first-pass analysis, the simulation programs 
used theoretical calculations that assumed the Galster 
parameterization for with different multiplicative fac- 
tors. We determined via simulation the optimal factor 
for each Q 2 that provided the best agreement between 
the simulated polarization ratios and the experimental 
asymmetry ratios. Next we fitted the current world data 
|H S [3 S B S E3 D3 and our first-pass acceptance 
and nuclear physics corrected results for G 1 ^ to a Gal- 
ster parameterization with two free parameters. Then 
we repeated the simulations using new theoretical calcu- 
lations that assumed this modified Galster parameteriza- 
tion for G%. As in the first-pass analysis, we determined 
via simulation the optimal factor that provided the best 
agreement between simulation and experiment. The dif- 
ferences between these analyses were negligible, and the 



TABLE I: Estimated systematic uncertainties in Ag/g [%]. 
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results from the two simulation programs agreed to bet- 
ter than two percent. 

The estimated values of the systematic uncertainties 
are listed in Tabled A significant advantage of our exper- 
imental technique is that the scale and systematic uncer- 
tainties are small; the analyzing power of the polarimeter 
cancels in the polarization ratio, and the beam polariza- 
tion, Pi, also cancels as it varied little during sequential 
measurements of the scattering asymmetries. We mea- 
sured the beam polarization with a M0ller polarimeter 
[l8| . and changes in Pl were typically on the order of 
one to two percent. The helicity of the beam was re- 
versed at a frequency of 30 Hz to eliminate instrumental 
asymmetries. 

A false asymmetry or a dilution of the asymmetry may 
arise from the two-step process 2 H(e, e'p) 1 H + Pb(p, n); 
the contamination from this process was assessed by run- 
ning with a liquid hydrogen target. The contamination 
levels are negligible (< 0.3%) for x = ±40° and ±90° 
at all of our Q 2 points, and for \ — 0°, the contami- 
nation levels are ~0.3% and ~3% at {Q 2 ) = 1.13 and 
1.45 (GeV/c) 2 , respectively; accordingly, we have not 
corrected our (Q 2 ) = 0.45 and 1.13 (GeV/c) 2 data for 
contamination from this two-step process. The net cor- 
rection obtained for the analysis of all of the data for 
(Q 2 ) = 1.45 (GeV/c) 2 [viz., for X = 0°, ±40°, and ±90°] 
amounted to 1.3±0.1%. In addition to charge- exchange 
reactions in the lead curtain, the flux of neutrons enter- 
ing the polarimeter may be depolarized as a result of nu- 
clear interactions in the lead curtain. Depolarization pro- 
cesses were simulated in GENGEN using a spin-dependent 
multiple-scattering algorithm employing quasifree scat- 
tering from a Fermi gas. The effects of depolarization 
cancel in the polarization ratio, and the residual non- 
cancellation effect upon g of less than 0.6% is included 
in the systematic uncertainty. 

Afanasev et al. calculated radiative corrections 
to the polarization-transfer coefficients, Ps'/Pl and 
Pl'/Pl- The primary effect is depolarization of the 
electron such that both polarization-transfer coefficients 
should be increased by ~1.9%, —3.7%, and —4.4% at 
(Q 2 ) = 0.45, 1.13, and 1.45 (GeV/c) 2 , respectively; how- 
ever, these effects nearly cancel in the polarization ratio 
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TABLE II: Results for g = G E /G n M and G%. 
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FIG. 3: The current world data on G E versus Q 2 ex- 
tracted from polarization measurements and an analysis of 
the deuteron quadrupole form factor 0, Q H H H HOl . The 
Galster parameterization §§ is the dashed line, and our two- 
parameter Galster fit (see text) is the solid line. 



such that the net effect upon g is small at (Q 2 ) — 0.45 
(GeV/c) 2 and negligible at the two higher Q 2 points. 

The values of g and G E that we report are listed in 
Table 1771 To determine our values for G E , we used the 
best- fit values for G 7 ^ (listed in Table [H} obtained using 
the methods described in [2(j • The new fit omits the data 
of [2ll |22| and includes the recent data of Xu et al. [2j| . 
The results are similar to those obtained by Friedrich 
and Walcher [24[ using a parameterization designed to 
investigate the role of the pion cloud. The fit is also 
similar to that of Kubon et al. [25| within their range of 
Q 2 , but extends to larger Q 2 . For g and G E , the first set 
of errors is statistical, and the second set is systematic. 
The quoted systematic uncertainties for g include a 2% 
uncertainty that results when a different selection of runs 
is used for the time calibration. To obtain the systematic 
uncertainties for G E , the relative uncertainties in G 7 ^ 
were added quadratically to the systematic uncertainties 
listed in TableUJand the 2% time calibration uncertainty. 

Our values for G E are plotted in Fig. together with 
the current world data on G E [J, Q, 0, la, la m 111 ex- 
tracted from polarization measurements and an analysis 
of the deuteron quadrupole form factor. We fitted these 
data and the G E slope at the origin as measured via low- 
energy neutron scattering from electrons in heavy atoms 
[l7| to a Galster parameterization: G% = — cl^tCd /{^-+ 
br), where r = Q 2 /AM 2 , G D = (1 + Q 2 /A 2 )- 2 , and 
A 2 = 0.71 (GeV/c) 2 . Our best-fit parameters are a = 
0.888 ± 0.023 and b = 3.21 ± 0.33. 




FIG. 4: Predictions of selected models (see text) for 
uM E /G p M and ^G^IGl, compared with proton [HEftl!! 
and neutron Q H H HH data. The neutron data 
symbols are the same as in Fig. 



Polarization measurements of G P E /G P M |2ll27ll2ll29l| 
and G E /G n M 0, S E 0, H are compared with pre- 
dictions of selected models in Fig. 0] The chiral soli- 
ton model |30j | reproduces the dramatic linear decrease 
observed in n P G p E /G p M for 1 < Q 2 < 6 (GeV/c) 2 ; how- 
ever, this model fails to reproduce the neutron data at 
large Q 2 . The light-cone diquark model [UJ achieves 
qualitative agreement with the low Q 2 proton and neu- 
tron data; however, at high Q 2 , it lies below(above) the 
proton(neutron) data. A calculation using the point- 
form spectator approximation (PFSA) with pointlikc 
constituent quarks and a Goldstone boson exchange in- 
teraction fitted to the meson and baryon spectrum |32l | 
also achieves qualitative agreement with the low Q 2 pro- 
ton and neutron data; however, it, too, fails to describe 
the high Q 2 proton and neutron data. A light-front calcu- 
lation using pointlike constituent quarks surrounded by 
a cloud of pions |33| describes the neutron data, but falls 
below the proton data at high Q 2 . A one-gluon exchange 
light-front calculation using constituent quark form fac- 
tors fitted to Q 2 < 1 (GeV/c) 2 data |34| agrees with the 
neutron data, but deviates from the proton data above 
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Q 2 ~ 3.5 (GeV/c) 2 . Finally, fits that couple vector me- 
son dominance with the predictions of perturbativc QCD 
|3fij agree with the entire range of the proton data, but 
fall below the neutron data above Q 2 ~ 1.2 (GeV/c) 2 . 

A successful model of confinement must be able to 
predict both neutron and proton electromagnetic form 
factors simultaneously. Although these models achieve 
qualitative agreement with the proton data, the neutron 
electric form factor is especially sensitive to small com- 
ponents of the nucleon wave function and differences be- 
tween model predictions for tend to increase rapidly 
with Q 2 . Our new G% data for 0.5 < Q 2 < 1.5 (GeV/c) 2 



provide a challenging test for confinement models and 
invite extensions to higher Q 2 . 
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